Shear flow and Kelvin-Helmholtz instability in superfluids 
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The first realization of instabilities in the shear flow between two superfluids is examined. The 
interface separating the A and B phases of superfluid ''He is magnetically stabilized. With uniform 
rotation we create a state with discontinuous tangential velocities at the interface, supported by the 
difference in quantized vorticity in the two phases. This state remains stable and nondissipative to 
high relative velocities, but finally undergoes an instability when an interfacial mode is excited and 
some vortices cross the phase boundary. The measured properties of the instability are consistent 
with a modified Kelvin-Helmholtz theory. 



Instabilities in the shear flow between two layers of 
fluids belong to a class of interfacial hydrodynamics 
which is attributed to many natural phenomena. Exam- 
ples are wave generation by wind blowing over water , 
the flapping of a sail or flag in the wind , and even 
flow in granular beds ||^ . In the hydrodynamics of invis- 
cid and incompressible fluids the transition from calm to 
wavy interfaces is known as the Kelvin-Helmholtz (KH) 
instability |^,^- Since Lord Kelvin's treatise in 1871, dif- 
ficulties have plagued its description in ordinary fluids, 
which are viscous and dissipative. They also display a 
shear-flow instability, but its correspondence with that 
in the ideal limit is not straightforward. The tangential 
velocity discontinuity in the shear-flow instability is cre- 
ated by a vortex sheet. In a viscous fluid a planar vortex 
sheet is not a stable equilibrium state and not a solution 
of the hydrodynamic equations 0. 

Superfluids provide a close variation of the ideal invis- 
cid limit considered by Lord Kelvin and thus an envi- 
ronment where the KH theory can be tested. The ini- 
tial state is a non-dissipative vortex sheet - the interface 
between two superfluids brought into a state of relative 
shear flow. So far the only experimentally accessible case 
where this can be studied in stationary conditions, is the 
interface between ^He-A and ^He-B where the order 
parameter changes symmetry and magnitude, but is con- 
tinuous on the scale of the superfluid coherence length 
^ ~ 10 nm. We discuss an experiment, where the two 
phases shde with respect to each other in a rotating cryo- 
stat: ■^He-A performs solid-body-Uke rotation while '^He- 
B is in the vortex-free state and thus stationary in the 
laboratory frame. While increasing the rotation velocity 
n, we record the events when the AB phase boundary be- 



comes unstable - when some circulation from the A-phase 
crosses the AB interface and vortex lines are introduced 
into the initially vortex-free B phase. On increasing the 
rotation further, the instability occurs repeatedly. Such 
a succession of instability events can be understood as a 
spin-up of ■^He-B by rotating ■^He-A. 

Our experimental setup is shown in Fig. 0. The 
AB boundary is forced against a magnetic barrier in a 
smooth-walled quartz container, by cooling the sample 
below Tab at constant pressure in a rotating refriger- 
ator. The number of vortices in both phases is inde- 
pendently determined from the simultaneously measured 
nuclear magnetic resonance (NMR) spectra |^,|l^. By 
tuning the barrier fleld or the temperature, the state of 
the sample can be changed from all A phase to all B 
phase or to a two-phase configuration. The evolution of 
the quantized vorticity as a function of Q is then observed 
to be radically different when the AB interface is present. 

The quasi-isotropic ^He-B supports singly quantized 
vortices with a core size comparable to f_j9|. In the 
anisotropic "^He-A we form vortex lines |10(| with con- 
tinuous skyrmion topology: Inside its central part the 
order-parameter amplitude remains constant but the axis 
of the orbital anisotropy covers a solid angle of Att. Such 
a "soft-core" structure carries continuous vorticity with 
two circulation quanta and is three orders of magnitude 
larger than the core of the B-phase vortex. Thus con- 
verting an A-phase vortex into a B-phase vortex requires 
large concentration of the flow energy. 

The large difference in core radii is the origin for the 
much lower rotation at which vortices start forming in A 
phase at the outer sample circumference, compared to B 
phase Pi . If the sample consists of only B phase then the 
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FIG. 1. Stabilization of the first-order ^He-A - ^He-B 
phase boundary. At pressures p > 21 bar, A-phase extends 
to lower temperatures in external magnetic field (see phase 
diagram on the right). The sample (length 11cm, radius 
R — 0.3 cm) is first cooled to A phase. On further cool- 
ing, the A^B transition happens in the coldest place at 
the bottom. The AB boundary then starts to move up and 
rises to a height z where the barrier field H{z) equals the 
value of the thermodynamic A— i-B transition Hab{T,p). The 
A phase in the top section remains in a metastable super- 
cooled state ||ll|. Ultimately, the boundary disappears when 
Hab{T,p) > [^(2)]max- For p — 29.0 bar and H{z) as shown 
in the figure with a current of 4 A in the barrier solenoid, the 
stable AB boundary exists below 2.07 mK down to 1.33 mK. 
The NMR spectrometers operate in homogeneous static mag- 
netic fields of 10 mT and 35 mT, chosen for best measuring 
sensitivity of the single-vortex signal. 



critical velocity for forming the first vortex line in the 
setup of Fig. |l| is VcB > 7mm/s |12 , while the critical 
velocity VcA is a factor of 20 smaller independently 
of the presence of the AB boundary. 

Vortex lines are thus easily created at low rotation in 
the A-phase section of the sample, while no vortices are 
detected in the B phase section. This is the ideal non- 
dissipative initial state where the two superfluid phases 
slide along each other without friction. When H. is in- 
creased further, sudden bursts of vortex lines are ob- 
served in the B-phase section, as shown in Fig. ^. The 
onset ftc depends on temperature and on the current in 
the barrier magnet. Overall, the measured characteristics 
of the bursts fit a shear flow instability, which provides a 
mechanism for the circulation to cross the AB interface. 

•^He-A and ■^He-B are states of the same order- 
parameter manifold. One of the conditions on their phase 
boundary is that the phase of the order parameter has 
to be continuous : If the circulation is not continuing 
across the interface ( Fig. ^A) then the existing A-phase 
vortex lines have to bend and form a vortex layer on the 
AB interface. This is the only hydrodynamically stable 
state, with vortex-free flow seen by the B-phase spec- 
trometer and a large cluster of vortex lines detected by 
the A-phase spectrometer. The build up of a sheet-like 
vortex layer means that the A-phase circulation cannot 
easily penetrate into the vortex-free B-phase and the AB 
interface remains stable up to the measured critical ve- 
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FIG. 2. Instability events during slowly increasing rota- 
tion {dfl/dt = 2.5 ■ 10^* rad/s^). NMR absorption signals 
as a function of are recorded (shifted arbitrarily on the 
vertical scale). The vertical jumps mark events in which 
SN new vortex lines enter the B phase section of the sam- 
ple. The height of each jump is proportional to SN, a small 
stochastic number. When averaged over a large interval of 
f2, the number of B-phase vortex lines grows linearly with Q. 
This is demonstrated also by the dashed line: the instabil- 
ity occurs independently of Q at constant critical drive, ie. 
Uc = |waB — Vn\r=R = ^IcR = coust. The three top most signal 
traces were recorded at the same temperature, ie. measure- 
ments at fixed T yield the same Sic (if small variations in 
VcA are accounted for [^). The two bottom traces at dif- 
ferent T illustrate that Qc depends on temperature. Here 
Tab{H = 0) = 0.785 Tc 



locity Uc = flcR ~ 2 - 4mm/s < UcB- 

At > r^c, some vortex lines have broken through the 
AB phase boundary. Thus there exists also a stable con- 
figuration in which vortices from the A phase continue 
into the B phase, where they form a cluster in the cen- 
ter (Fig. 1^). The likely topology of the intersection is 
illustrated in Fig. ^C. Thus, although the AB interface is 
not directly monitored by the two NMR spectrometers, 
the hydrodynamically stable states below and above flc 
can only be the configurations in Fig. ^. By comparing 
B-phase vortex-line creation in the presence and absence 
of the AB interface, we conclude that the events in Fig. ^ 
originate from the AB phase boundary. 

In the classical KH instability the interface between 
fluids with densities pi and p2 becomes destabilized by 
inertial effects, which are normally balanced by gravity g 
and surface tension a. The relative velocity |w2 — | acts 
as a drive. When it reaches a critical value 1^ given by 



PlP2 

Pi + P2 



■{V2 - ViY 



(1) 



where F = g{pi — P2) is the gravitational restoring force, 
waves with wave vector k = F/a are created on the 
interface. 

This approach can be generalized for superfluids in 
terms of two- fluid hydrodynamics |15|. For superfluid 



2 



n>Qc 




FIG. 3. The two configurations of vortex fines with the 
AB boundary in rotation. (A) Q < fl^: Vortex lines are 
formed in '^He-A while '^He-B remains vortex-free. Near the 
AB boundary the A-phase vortices bend paraUel to the in- 
terface and form a vortex sheet between the sliding super- 
fluids. The radial distributions of the normal and superfiuid 
velocities far from the AB interface are depicted below, with 
VnA ~ VnB ~ , while v^B ~ 0. (B) Q > Q^'- Vortices are 
observed to appear in the B phase in events of a few lines 
at a time. They form a central cluster in the B-phase sec- 
tion. (C) A hydrodynamically stable state with respect to 
externally imposed perturbations in fi, T, or H exists at the 
AB interface for > i^c. A topologically stable configura- 
tion for the vortex-line intersection with the AB interface is 
suggested in Ref. |l4|: The doubly quantized A-phase vor- 
tex terminates at the AB interface in two point singularities, 
known as boojums. These, in turn, are the end points of two 
singly quantized B-phase vortices. 

•^He under our experimental conditions it is safe to as- 
sume that the normal fractions are always in solid-body 
rotation. Instead of gravity, the restoring force is now 
produced by the magnetic barrier H{z), owing to the 
difference in the susceptibilities xa and xb {T, H) : F — 
(1/2)(XA - XB(r,i/))V(iJ2) at HABiT). The restoring 
force can thus be calculated as a function of temperature 
and current in the barrier solenoid. The motion of the 
AB interface with respect to the normal component is 
subject to a finite damping p6| . Since the initial state 
is non-dissipative, damping does not explicitly appear in 
the stability condition, but modifies the build-up rate of 
the interface perturbations when the instability develops. 
The onset of the instability can be derived from the dy- 
namics of small amplitude perturbations, as Eq. (|l|) is 
generally derived in textbooks |Q, or from the thermo- 
dynamics when perturbations of the interface lead to a 
negative free energy in the rotating frame. The mode for 
which the interface first becomes unstable has the same 
wave vector k = ^ F ja as before at a drive given by 

^PsA (WsA - Wn)^ + ^PsB (t^sB " W„)^ = \faF , (2) 

where Ps, Pm and Wgj are the densities and velocities of 




FIG. 4. Critical velocity fic for the first appearance of 
B-phase vortex lines in the presence of the AB boundary 
as a function of temperature, while is slowly increased 
(f2 = 5 ■ 10~^rad/s^). The current in the barrier solenoid 
is constant but different for the three sets of data. It con- 
trols the magnetic force F which is largely responsible for the 
shape of the curves. The solid curves represent Eq. (^) if one 
sets f2c = |«sB — Vri\IK and WsA = Vn- No fitting parameters 
are used. The values for cr{T), xa ~ XbCF H), and Pe{T) are 
obtained from accepted references ^,|l8|-po| , while H and S7H 
apply for the profile H{z) at the location of the AB boundary: 
Hab{T)=H{z). 



superfiuid and normal components [0. In Fig. ^ we plot 
the prediction for the instability for three magnetic bar- 
rier profiles. With no fitting parameters the agreement 
with experiment is surprisingly good. 

From Fig. || it is seen that of order 6NRi 10 vortex lines 
enter the B phase in one instability: in fact, at 0.77 Tc 
(with Uc — 0.39 cm/s at 29.0 bar) the measured average 
for more than 100 events is 6N 11. To compare with 
Eq. (||), SN corresponds to the number of circulation 
quanta k = h/2m^^ which fit in one corrugation of the 
interface mode of size Ac/2 = 7r/fc. In practice, the coun- 
terflow |wsA— ■^nl in A-phase is limited by the small critical 
velocity VcA (l^ , and we may set Vsa — f n ~ in Eq. (||) . 
In solid-body rotation there are then N ruTiR^ Vic/ n vor- 
tex lines in the A phase which all flare out into the lateral 
sample boundary at the AB interface (Fig. |3|A). Mea- 
sured along the perimeter of the sample there are Oci?//c 
circulation quanta per unit length and thus in one cor- 
rugation 6N K {ttUc) / {kn). From Eq. (^) this is seen to 
be SN Ri {27ra)/ (kUcPsb), which at 0.77rc gives 9 vortex 
lines in agreement with the measured number. 

Eq. (H) has more interesting properties. As the normal 
components of both phases are in solid-body rotation and 
thus at rest in the rotating frame, it is the "superfluid 
winds" - the flow of the superfiuid component, \vs — Wn|, 
on each side of the interface - which produce the insta- 
bility. It takes place even if the two superfluids have 
the same densities and velocities. In this sense it resem- 
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bles the flapping flag instability discussed by Rayleigh 
1^: The role of the flagpole, which fixes the reference 
frame, falls here on the normal component, which moves 
with the rotating container. 

It is interesting to compare Eqs. (|l]) and (||) in the low 
temperature limit, when pn ^ and psA ~ PsB ^ P- We 
find that Eq. gives a critical velocity which is by \/2 
smaller than that from Eq. In fact, at T = the 
valid condition is Eq. since pn is exactly zero and no 
normal component is left to provide a reference frame for 
the now undamped interface The ideal KH condi- 

tion is then restored, as the only remaining drive is the 
velocity difference between the supcrfluid components. 
Not surprisingly, the two cases Pn — and p^ ^ give 
different onsets for the instability: For a flnite-size sys- 
tem the result would depend on the observation time - 
the time one waits for the interface to be coupled to the 
laboratory frame. 

The shear-flow instability in superfluids thus provides 
an explicit bridge between the ideal inviscid and the vis- 
cous hydrodynamics. For two-fluid hydrodynamics it 
may have important consequences. It represents a new 
surface mechanism for vortex generation which might be 
applicable also to other interfaces, of which the normal 
state - superfluid boundary is of major interest. This in- 
terface occurs, for example, when the bulk superflow in- 
stability limit is reached in rotating '^He-B at the sharpest 
surface spike on the lateral sample boundary |^ . A bub- 
ble of normal liquid is then created around this spike. 
Since the surrounding superflow moves at a high relative 
velocity, the bubble's interface may undergo a shear-flow 
instability and promote the generation of a vortex line. 

A related case has been investigated in bulk ^He-B 
which is in vortex-free rotation while it is irradiated with 
slow neutrons |pl|| . The absorption process of one neu- 
tron creates locally a bubble of 100 fim size which is 
overheated up to the normal state. The bubble then 
quench cools back to the temperature of the surround- 
ing superflow, as the normal-superfluid interface of the 
bubble contracts at high velocity. The measurements 
prove directly that vortices are formed in this rapid non- 
equilibrium transition and that the process is consistent 
with a volume effect known as the cosmological Kibble- 
Zurek (KZ) mechanism p^ , p3| . It has been suggested as 
the origin for cosmic-string formation when rapid non- 
equilibrium phase transitions swept across the Early Uni- 
verse. The ^He-B neutron experiment has been exam- 
ined by solving analytically and numerically the time- 
dependent Ginzburg-Landau equation p^ . This study 
shows that in the non-equilibrium transition of the neu- 
tron bubble there exists also a competing surface mecha- 
nism, when the rapidly moving normal-superfluid inter- 
face becomes unstable and vortex rings form around the 
bubble in a corrugation instability. This is direct theoret- 
ical evidence for the existence of the shear-flow instabil- 
ity also at the normal-superfluid interface, between the 



surrounding superflow and the stationary normal liquid 
within the bubble. A similar situation is the normal- 
superfluid boundary in a gaseous rotating Bose-Einstein 
condensate. Here the shear-flow instability also leads to 
a deformation of the surface of the condensate cloud by 
ripplons and to vortex-line formation [ p5| . 

The shear-flow instability thus appears to be a univer- 
sal surface mechanism for defect formation in superfluids 
and as intrinsic as the KZ mechanism in the bulk liq- 
uid. We are left to wonder whether or not these two are 
the only possible processes in the clean limit at the low- 
est temperatures, to which other observations of vortex 
formation can be reduced. 
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